
#include "Interpolation.hpp"

int main(){
    std::cout<<"The Assignment D tests Hermite Interpolation."<<std::endl<<std::endl;
    std::vector<double> x{0,0,3,3,5,5,8,8,13,13},y{0,75,255,77,383,80,623,74,993,72},d_o{0,1,0,1,0,1,0,1,0,1};
   Hermite p(x,y,d_o);
    std::cout<<"Hermite Interpolation p(x)="<<p.getPol()<<std::endl;
    std::cout<<"Hermite Interpolation p'(x)="<<p.getPol().diff()<<std::endl;
    std::cout<<"(a) When t=10s, p(t)=p(10)= "<<p.getPol()(10)<<std::endl;
    std::cout<<std::endl<<"--------------------Test succeeds.--------------------"<<std::endl<<std::endl;
    return 0;
}

